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Abstract 

In a recent paper, we have reported a universal power law for both 
site and bond percolation thresholds in any Bravais lattice with q equiv- 
alent nearest neighbors in dimension d. We now extend it to three differ- 
ent classes of lattices which are, respectively, anisotropic lattices whith 
not equivalent nearest neighbors, non-Bravais lattices with two atom 
unit cells, and quasicrystals. The investigation is focussed on d = 2 and 
d = 3, due to the lack of experimental data at higher dimensions. The 
extension to these lattices requires the substitution of q by an effective 
(non integer) value q e ff in the universal law. For each out of 17 lattices 
which constitute our sample, we argue the existence of one q e ff which 
reproduces both the site and the percolation threshold, with a deviation 
with respect to numerical estimates which does not exceed =f0.01. 



*Laboratoire de FUniversite P. et M. Curie - Paris 6, associe au CNRS (URA n° 800) 



1 



1 Introduction 



Very recently we have posulated a universal power law for both site and bond 
percolation thresholds [1]. The formula yields thresholds for any Bravais lat- 
tice, at any dimension with an impressive accuracy. It writes 

Pc = Po [(d-l)(q-l)]- a d b (1) 

with d the space dimension and q the coordination number. While b = a 
for bond dilution, it is b = for site dilution. Three different classes were 
found with three distinct parameter sets {p ; a}. The first class includes two- 
dimensional triangle, square and honeycomb lattices. It is characterized by 
{po = 0.8889; a = 0.3601} for site dilution and by {p = 0.6558; a = 0.6897} 
for bond dilution. Two-dimensional Kagome and all other lattices of cubic 
symetry (for 3 < d < 6) constitute the second class which is characterized 
by {p = 1.2868; a = 0.6160} and {p = 0.7541; a = 0.9346} for sites and 
bonds respectively. A third class has been found at high dimensions (d > 6), 
which recovers the infinite Cayley tree limit, but is not relevant to present 
investigation which deals with lattices only at d — 2 and d = 3. 

Besides the dimension d, percolation thresholds within a class depend 
only on q. This is understood in lattices where the q nearest neighbors of any 
site are equivalent, which is indeed the case for all lattices above mentioned. 
This is however a drastic restriction since many percolation problems in physics 
deal with lattices which do not have this property. It is the purpose of the 
present work to investigate the extension of Eq. (1) to other lattices, via the 
substitution of q by an effective parameter q e ff- 

We have checked on a sample of 17 lattices that, for each of them, there 
exists one value of q e ff which reproduces both the site and the bond percolation 
thresholds. The error is within ±0.006 for all the lattices, except for two 
of them (dual of Penrose, and dodecagonal lattice with ferromagnetic links), 
where the error reaches 0.01 for reasons we also discuss in this paper. We 
consider this sample is large enough to be representative of a general trend, 
and then conclude in the existence of such a q e ff in any lattice. For a lattice 
which does not belong to our sample, this parameter could be used as an 
intermediate quantity to predict bond from site percolation thresholds, or vice 
versa, with the same accuracy. 
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The paper is organised as follows. Section 2 discusses a series of lattices 
with non-equivalent neigbhors. The introduction of an effective number of 
nearest neighbors in our universal formula for percolation is argued in Section 
3. The results are shown and discussed in Section 3. The last section contains 
some remarks for future work. 

2 Non-equivalent neigbhors 

Some periodic Bravais lattices are anisotropic, for instance the hexagonal lat- 
tice at d — 3. In this case, any lattice site has 6 equivalent nearest neighbors in 
the a, b plane (bonding angle is 60°) and 2 non-equivalent sites along the c axis 
(bonding angle is 90°). Actually, the anisotropic lattice percolation threshold 
should depend on the degree of anisotropy. For the hexagonal lattice, it means 
p c should be different from that of the corresponding isotropic lattice with the 
same set d — 3, q — 8, i.e., the bcc lattice. This difference has been observed 
indeed recently [2] in the particular case of the stacked triangular lattice, which 
becomes the hexagonal lattice when a = b = c. 

There also exist non-Bravais lattices which are periodic like fee and 
bcc. Another case is the hexagonal close packed (hep) lattice which, on a 
topologic view point, is a simple hexagonal lattice with two atoms per unit 
cell. Percolation thresholds for the hep were obtained long ago [3]. 

Some lattices are not even periodic. This is the case of quasicrystals 
which are aperiodic lattices with long-range order. Besides their own interest 
such structures can serve as models for alloy-like materials making a growing 
interest in the determination of the quasicrystalline lattice percolation thresh- 
olds. First determinations of percolation thresholds on 2-d quasilattices have 
been made on the Penrose tiling and its dual [4, 5, 6]. Recently, percolation 
thresholds have also been computed in two of the most important quasilattices, 
the simple octagonal and dodecagonal tilings [7], which belong to the Penrose 
local isomorphism class. 

By analogy with other lattices, one can require that lattice sites are 
connected only via the tile edges. The corresponding percolation problem is 
named chemical percolation. Adding connection through the diagonals of the 
tiles which are shorter than the tile edge leads to the so-called ferromagnetic 
percolation [8]. Both these percolations will be considered here. These lattices 
are reproduced in Fig. (1). 

At last, some lattices do not have a single- valued coordination number q. 
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Lattices with mixed- valued coordination can be either periodic or non-periodic. 
An exemple of periodic lattice at d = 2 is provided by the dice lattice which 
mixes q = 3 and 6. An exemple of non-periodic lattice is the Penrose tiling, 
which mixes q = 3, 4, 5, 6 and 7. Note the dual of dice is the Kagome lattice 
and vice versa, while the dual of the Penrose tiling is named dual lattice of 
Penrose. Both the Kagome and the dual of Penrose are lattices with the 
single-valued coordination q — 4. Octagonal and dodecagonal lattices have 
mixed-value coordinations. Numerical estimates of the percolation thresholds 
of dice can be found in [6]. 

3 An effective number of nearest neigbhors 

In the past, attempts have been made to generalize empirical relations such 
as the Scher and Zallen approximate [9] which depends on dimension. For 
instance at d — 2, these are qp c (bond) ~ 2.0 and fp c (site) ~ 0.45, with / 
the lattice filling factor. Extension to quasilattices [6, 7] consists in replacing 
q by mean coordination number z in the invariant appropriate to bonds (the 
problem for sites is not solved, because there is no ambiguous definition of the 
filling factor in quasicrystals [7]). 

In the same spirit, we propose to extend Eq. (1), replacing q by an 
effective value q e ff- The dependence of the percolation thresholds on site con- 
nectivity does not imply that the relevant variable should be the arithmetic 
average z. Therefore, we regard q e ff as a parameter which has to be deter- 
mined, rather than arbitrarily forced to be equal to z. 

The need of substituting q by q e ff in the equations is easily understood 
in both the case of quasilattices and the case of lattices which mix different 
values of q. It is also needed in anisotropic lattices. Let us again consider the 
hexagonal lattice with spacing a, b and c in the three lattice directions. In 
the limit c — > +oo, one is left with decoupled ab planes for physical systems 
having a finite range of interaction. Therefore, the percolation threshold of 
the hexagonal lattice must depend somehow on the ratio c/a. It is expected 
to range between that of isotropic bcc with d = 3, q = 8 and that of the d = 2 
triangular lattice, associated to the limit c/a — > +oo. Note in this limit one 
recovers an isotropic lattice with q = 6 instead of q — 8, which suggests an 
effective coordinance 6 < q e ff < 8. In view of such considerations, we propose 
to extend Eq. (1) substituting q by some q e ff for any lattice. 

Note within the two classes defined by the set (a,po), q e ff is the unique 



4 



unknown parameter. For each lattice, we find a value of this fitting param- 
eter which, when inserted into Eq.l, reproduces both the site and the bond 
percolation threshold within =(=0.01. Values of q e ff are reported for different 
lattices in Table (1), together with site and bond percolation thresholds p s c , 
p\ obtained when q is replaced by q e ff in Eq. (1). Mean coordination z and 
pi for site and bond percolation are also reported for comparison. Exact or 
numerical estimates p e c are from Refs [6, 7, 10]. 

4 The results 

For lattices at d — 2, q e ff differs from z by 1% in the case of the dual latttice 
of Penrose, and is smaller than 0.5% for all the other lattices. These data 
corroborate that two dimensional lattices are divided into the two distinct 
classes defined in [1]. This is illustrated in Figs. (2) and (3), where p^s are 
reported in a log — log plot such that the experimental points for lattices in 
the same class align on a straight line according to Eq.l. 

The pertinent variable after Eq. (1) is (d — l)(q e ff — 1) for sites and 
(d — l)(q e ff — l)/d for bonds. Those are the variables in abscissa in Fig. (3) 
which reports data for lattices in dimensions d = 2 and d = 3 (data in higher 
dimensions have been already displayed in a similar plot in [1] with q e ff = q). 
Since the first class concerns only lattices which are all at d — 2, above variables 
for respectively sites and bonds may reduce to one single common variable q e ff 
as shown in Fig. (2). | A | reaches 0.01 only in the dodecagonal lattice with 
ferromagnetic links, and in the dual lattice of Penrose. In all the other cases, 
the error in the percolation threshold estimates is only on the third decimal. 

Note the larger error in the dodecagonal lattice. We attribute it to an 
actual bond percolation threshold p e c {bond) = 0.495 [6] larger than a priori 
expected. One would indeed have expected p c (bond) = 0.475 from the straight 
line in Figure 3, in which case a value q e ff = 4.289, close to z — 4.27, would 
have reproduced both site and bond percolation thresholds within | A | = 
0.001. The problem with the estimate p e c (bond) = 0.495 is also evident from 
the Scher and Zallen invariant. It yields zp e c (bond) = 2.11, the largest value 
of this parameter among all the lattices investigated [6]. At the opposite, 
p c (bond) = 0.475 would yield zp c (bond) = 2.028, close to the invariant value 
2.0 at d = 2. 

At d — 3, percolation values are reproduced with a very good accuracy, 
since | A |< 0.004 for all the lattices. The difference between q e ff and z is not 
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negligible as it reaches few per cent in some cases. This is easily understood in 
the case of the hexagonal (stack triangle) case, where we have argued earlier 
that a value q e ff smaller than 8 is expected. Actually, we find q e ff = 7.66, 
as a consequence of the anisotropy. The hexagonal close packed (hep) lattice 
has percolation thresholds which are different from those of the fee lattice, as 
expected since the hep lattice is not a Bravais lattice. 

However, differences are small which may be attributed to the fact that 
both lattices are indeed isotropic, each site being surrounded by 12 equivalent 
neighbors. In this context, the small value of q e ff close to 11 in hep lattice is 
not only due to the non-Bravais nature of the lattice. It is also related to the 
fact that in the fee lattice, q e ff is only 11.6, significantly smaller than q = 12, 
although the coordination number is single-valued. 

On the other hand, q e ff differs from q by only 1% in the other lattices 
(sc, diamond, bee). Yet one would expect q e ff = q in such isotropic lattices 
with single- valued coordination number. The difference between q e ff and q in 
this case illustrates that our formula for the percolation thresholds is not exact 
as we already stated in [1], and shown convicingly in ref. 2. Nevertheless, both 
site and bond percolation thresholds for all the lattices in any dimension are 
provided within 1% by our universal law involving only two parameters: the 
dimension d and a parameter q e ff which contains more information on the 
geometry than the mean coordination. 

5 Conclusion 

We have shown that our universal formula for percolation threholds we re- 
ported earlier [1] for periodic Bravais lattices with equivalent nearest neigh- 
bors does extend to any kind of lattice, provided the coordination number is 
replaced by an effective value g e //. We then conclude that a good estimate of 
both site and bond percolation thresholds can be obtained from the formula 
in Eq. 1 involving only the dimension d of the lattice, and one parameter 
q e ff which contains the geometric information of the lattice. This parameter 
however, does not reduce to the mean coordination number z, although z and 
q e ff differ by few per cent only. Indeed, we find the universality does include 
different numerical values for q e ff and thus different percolation thresholds for 
lattices which have the same set (d, z). This is evident from Table (1) which 
reports results for as many as six lattices with d — 2, z — 4, two lattices with 
(d — 3, z = 8), two other ones with (d — 3, z — 12). 
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In our previous work, we had only three lattices belonging to the first 
class. The present work extends this class from 3 to 8 lattices. We might have 
invoked chance for three points aligned on a same line, but not for 8 points 
like in Fig.l. Therefore, we confirm the existence of two different classes, 
one for some of the two-dimensional lattices, the other one for all the other 
two-dimensional lattices and all the lattices up to d=7. 

We do not have a scheme to derive the relevant variable q e ff from the 
geometry of the lattice. However, an important result of this work is that this 
variable does exist. It means that there exists a single value for q e ff which 
accounts for both the site and the bond percolation thresholds for any given 
lattice within =(=0.01. This result is sufficient to give our formula a prediction 
ability for pecolation thresholds of other lattices which have not been computed 
yet. For example, the knowledge of one site (bond) percolation threshold for 
a given lattice is sufficient to determine a point on the relevant straight line in 
Fig. 2 or 3. Then q e ff can be found from the absissa of this point, which in 
turn allows for the determination of the bond (site) percolation threshold from 
our universal formula, within =(=0.01. Depending on whether p s c or p\ is known, 
the value of q e ff deduced from Fig. (2) or (3) will be different. However, 
this difference only corresponds to the deviation of the p^s with repect to the 
universal law, i.e. less than 1% with the exception of only few out-liers. 

The robustness of our formula suggests the extension to more complex 
problems such as directed percolation. Also, anisotropic percolations with 
different bond probabilities in different directions may be adressed in the near 
future. 
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Figure captions 



Fig. 1. Less common lattices studied in this paper : Kagome (a) and 
its dual, Dice lattice (b); Penrose quasicrystalline lattice (c) and its dual (d); 
octagonal quasicrystalline lattice with ferromagnetic links (e) and chemical 
links (f); dodecagonal quasicrystalline lattice with ferromagnetic links (g) and 
chemical links (h). 

Fig. 2. Inverse of percolation thresholds as a function of the variables 
(d — l)(q e ff — 1) and (d — l)(q e ff — l)/d) which reduce here (d = 2) to the 
single variable q e ff in logarithmic scales for lattices belonging to the first class. 

Fig. 3. Inverse of percolation thresholds as a function of the variables (d— 
l)(q e ff — 1) and (d— l)(q e ff — l)/cf) approriate to site and bonds, respectively, 
in logarthmic scales for lattices belonging to the second class. 
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0.5928 
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+0.0042 
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0.4935 


-0.0064 


Honeycomb 
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2.99 
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-0.0024 


0.6527 
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+0.0054 


Triangular 
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0.5 


0.4986 


-0.0014 


0.34729 


0.34955 


+0.0023 


Dice* 
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4.189 


0.5851 


0.5854 


+0.0003 


0.476 


0.4754 


-0.0006 


Penrose* 


4 


4.194 


0.5837 


0.5851 


+0.0014 


0.477 


0.4748 


-0.0022 


Octa.chem. links* 


4 


4.170 


0.585 


0.5867 


+0.0017 


0.48 


0.4773 


-0.0027 


ct a . f erro . links * 


5.17 


5.013 


0.543 


0.5389 


-0.0041 


0.40 


0.406 


+0.0057 


Dode.chem.links* 


3.63 


3.638 


0.628 


0.6269 


-0.0011 


0.54 


0.5419 


+0.0019 


second class 
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Lattice 


z 
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Pc 
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P 

Pc 


Pc 


A 


Kagome 


4 


3.980 


0.6527 


0.6567 


+0.0040 


0.5244 


0.5195 


-0.0049 


dual Penrose 


4 


4.04 


0.6381 


0.6487 


+0.0106 


0.5233 


0.5099 


-0.01341 


Dode.ferro. links* 


4.27 


4.218 


0.617 


0.6264 


+0.0094 


0.495 


0.4835 


- 0.0115 


hex. compact 


12 


11.146 


0.204 


0.2015 


-0.0025 


0.124 


0.1263 


+0.0023 


stac. triangle 


8 


7.661 


0.2623 


0.2611 


-0.0012 


0.1859 


0.1872 


0.0013 


Diamond 


4 


4.0087 


0.43 


0.4260 


-0.0040 


0.1859 


0.3935 


+0.0055 


sc 


6 


5.9558 


0.3116 


0.3132 


+0.0016 


0.2488 


0.2468 


-0.0020 


bcc 


8 


8.1355 


0.246 


0.2502 


+0.0042 


0.1803 


0.1755 


-0.0047 


fee 


12 


11.626 


0.198 


0.1958 


-0.0022 


0.119 


0.1210 


0.0020 



Table 1: percolation thresholds from this work p c compared to "exact estimates" 
pi taken from [6, 7 and 10]. A = p c — p e c . * refers to multi-valued coordination 
number (mean coordination z). All the lattices are in dimension d — 2, except 
the 6 last ones in the second class which have d = 3. We have completed the 
table caption as follows: q e ff has been chosen as the arithmetic average of the 
parameters q which reproduces pi for respectively site and bond when inserted in 
Eq. (1). 
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